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Emergence of a Barchan Belt in a Unidirectional Flow: Experiment and Numerical 

Simulation 
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^ Department of physics, Osaka University, Osaka 560-0043 
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We observed time evolution of dune fields in a water tank experiment and simulated it by 
using a simple model without taking complex fluid dynamics into account. The initial sand bed 
changed its form into transverse ripples , that is, dunes with straight crest lines perpendicular to 
the flow direction. Then the crescent ic shaped dunes called barchans emerged from transverse 
ripples. 

KEYWORDS: sand dune, barchan, transverse ripple, dune field, water tank experiment, numerical simula- 
tion 



Sand dunes, which are found in deserts, on the sea 
bottom and also on Mars, have various morphologic pat- 
terns. Dunes are formed by interactions between the flow 
(wind or water) and sand. The flow makes the shape of 
a dune by transporting sand particles. The dune topog- 
raphy, in turn, acts as a boundary condition on the flow. 
Types of sand dunes are determined by the amount of 
available sand and the variation of the flow.^ When the 
flow is unidirectional and the amount of sand is abundant 
for covering the entire bedrock, dunes with straight crest 
lines perpendicular to the flow called transverse dunes 
are generated. When the flow is unidirectional but the 
amount of sand is insufficient, on the other hand, crescen- 
tic shaped dunes called barchans are formed. Barchans 
usually migrate as a group and form a barchan belt?'^ 
The barchan belt is generated from a sand source which 
exists windward.^ When there are transverse dunes wind- 
ward which act as a sand source, a morphologic sequence 
of three types of dunes, transverse dunes, barchanoids 
and barchans from upwind to downwind is observed in 
deserts, where barchanoid is a dune with a wavy crest. ^ 

Observing the whole process of the formation dynam- 
ics of a barchan belt is quite difficult in general because 
time scale of the dynamics is very long. To overcome 
this difficulty, laboratory experiments using a water tank 
were proposed, , which reduce the time scale of dy- 
namics by less than a hundredth in rough estimates. 
Quite recently, a transition from transverse dunes to 
barchans was observed in a water tank under an oscilla- 
tory flow by Endo et. alJ A water tank experiment under 
an unidirectional flow was also conducted and a similar 
morphologic transition was observed.^ 

Only a few studies on the barchan belt have been made 
theoretically or by numerical simulations. Lima et. al}^ 
and Hersen et. al.^ has studied the stability of a barchan 
belt, assuming that many barchans somehow formed. 
However, formation processes of the barchan belt from 
a sand source have not been discussed. Moreover their 
models have not taken into account a solitary wave be- 
havior of barchans, which was found in collision dynam- 
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Fig. 1. Height dependence of the flow velocity in a water tank 
experiment. 



ics.^^^^ In this letter, we will investigate formation dy- 
namics of a barchan belt by using a numerical model, 
which has succeeded in reproducing a solitary wave be- 
havior. We also conduct a water tank experiment for 
the purpose of observing development of a barchan belt. 

The experiment was performed in a water tank (10 m 
long, 20 cm wide and 50 cm deep) in which a pump can 
generate a steady unidirectional flow. Water depth 13.5 
cm was kept constant. The sand was well-sorted: 80-100 
fim in diameter and 130 g in total weight. In the initial 
condition, the sand bed was placed in an area of 15 cm 
long and 20 cm wide. Figure shows the flow velocity 
measured at several different heights in a run without a 
sand bed. Since height of the sand bed was less than 1 
cm, the transportation of sand was considered to occur 
inside the boundary layer. 

After the flow started, stages of morphologic transi- 
tion were observed. First, instability on the sand bed 
grew and transverse ripples started to appear (Fig. 
2(b)). Next, parts of the transverse ripples were breached 
at some places and the fragments eventually became 
barchans (Fig. 2(c) and 2(d) ). Figure 2(d) closely resem- 
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Fig. 2. Time evolution of a dune field in the water tank exper- 
iments under a steady unidirectional flow, (a) second (initial 
conflguration) , (b) 180 seconds, (c) 480 seconds and (d)1140 sec- 
onds. 

bles a barchan belt in a real desert sketched by Bagnold.^ 
Numerical simulation of a dune field taking into ac- 
count complex process of flow dynamics and granular 
dynamics is a formidable task because of computational 
cost. In what follows, we used a simple model with min- 
imal processes. Both complex fluid dynamical processes 
such as convection and viscosity and effect of size and 
shape of sand grain were ignored. This model has already 
succeeded in reproducing morphologic features of a single 
barchan and collision dynamics of two barchans.^^ 

In the model, the dune field is divided into square 
cells. -^^'-^^ Each cell is considered to represent an area 
of the ground which is sufficiently larger than the sand 
grains. A field variable h{x,y,t) which expresses the lo- 
cal surface height is assigned to each cell; t denotes the 
discrete time step and the spatial coordinate x and y are 
positions of the center of a cell in the fiow direction and 
in the lateral direction, respectively. The edge length of 
the cell is taken as a unit of length. In short, x, y and t 
are discrete variables while h{x,y,t) takes a continuous 
value. 

Next, we models a motion of sand grain considering 
only the following two processes: saltation and avalanche. 
Saltation is the transportation process of sand grains by 
the fiow. The schematic diagram of the saltation is shown 
in Fig. 3. The saltation length and saltation mass are de- 
noted Ls and respectively. Here, the saltation mass in- 
deed is a volume of sand transfered from a cell to another 
cell. Since the area of a cell is unity, q represents change 
of height of a cell by saltation. In each time step of a 
simulation, the saltation mass q shifts from a cell (x, y) 
to the leeward cell {x -\- Ls^y). Hence changes of height 
h{x, y) h{x, y)-q and h{x -\- Ls,y) ^ h{x-\- Ls,y)-\-q 
take place at the taking-off cell and the landing cell, re- 
spectively. Saltation length is then modeled by the fol- 
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Fig. 3. Schematic diagram of saltation process. Ls and q are 
saltation length and saltation mass, respectively. 
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Fig. 4. Relationship between the saltation length and the dune 
height (eq.(l)). Saltation length is rounded up to an integer. 



lowing equation: 

Ls = a + hh{x,y,t) - ch^{x,y,t). (1) 

The phenomenological parameters a=1.0, &=1.0 and 
c=0.01 are used in this work. Ls is rounded up to an 
integer as is shown in Fig. 4 . Equation (1) represents 
that the sand is transported farther away as height is 
higher, but not too far. Note that equation (1) is used 
only in the range in which Ls is an increasing function 
of h{x^y^t). The saltation mass q is fixed as 0.1 for sim- 
plicity throughout the present work. 

Avalanche, on the other hand, is the process in which 
the sand slides down along the steepest slope if it be- 
comes steeper than the angle of repose. In this simula- 
tion, first, the cells are marked if the slope, which makes 
with the lowest of their nearest-neighbor cells, exceeds 
the angle of repose. Then half of the excess volume of 
each marked cell is transfered. Here the angle of repose 
is fixed as 34°. This procedure is repeated until all cells 
satisfy the stability condition. Two processes of saltation 
and avalanche are performed by turns. 
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Fig. 5. Time evolution of the dune field for Do = 3 and Lq = 256 
in numerical simulation. The center of mass coordinate is fixed 
in the figure. Time steps are (b) 1000, (c) 2000 and (d) 5000. 



In the numerical field, the open boundary condition 
was imposed in the flow direction while the lateral 
boundary was impenetrable to sand grains. It means that 
there was no influx sand from the upwind boundary and 
sand which went out from the downwind boundary was 
simply deleted. The dune field consisted of 1052 (wind- 
ward) by 512 (lateral) cells. 

As an initial condition, we set a rectangular sand bed 
in order to reproduce the experimental situation. The 
initial sand bed had a sand depth h{x^ 0) = Dq = 3, 
length Lq = 256 along the flow direction and width l^o = 
512 (Fig. 5(a)). Height of the sand h{x^y^t) were then 
perturbed with ± 10 % random noise. 

After the simulation started, instability grew in the 
flow direction and made the sand bed change to trans- 
verse ripples. After a while, barchans emerged from the 
transverse ripples because of instability in the lateral di- 
rection (Fig. 5) . The processes quite similar to what 
were observed in the water tank experiment were thus 
reproduced. 

Next, we changed the sand depth I^o and the length 




Fig. 6. Dependence of ripple size on the depth Dq. The length 
and time steps were fixed as Lq = 512 and t = 10000, respec- 
tively, (a) Do = 2, (b) Do =4 and (c) Do = 6. 



Lq of the initial sand bed. As I^o became shallower, rip- 
ples were found to be narrower (Fig. 6). As Lq became 
larger, on the other hand, the width of ripples was un- 
changed but the number of ripples increased (Fig. 7). It 
means that the ripples have a characteristic length which 
depends on initial depth of the sand bed as long as I^o 
is much larger than q. 

We observed two types of transition from transverse 
ripples to barchans. Type A, an example of which is 
shown in Fig. 8, occured when a portion of a transverse 
ripple happened to be much lower than the rest body. 
The transverse ripple gradually extended to the leeward 
at that part and eventually breached because the migra- 
tion velociy is inversely proportional to the height. ^^^^ 
After that, the fragment thus formed became a barchan. 
Transition of this type were also observed in the water 
tank experiment.^ Type B, an exapmle of which is shown 
in Fig. 9, occured when a portion of a transverse ripple 
happened to be much higher than the rest body. The 
transverse ripple formed a crescent shape at that part, 
which evetually breached from a ripple to form an iso- 
lated barchan. In this case, a barchan appeared in the 
windward of the ripple. This type was observed only by 
the simulation so far. Whether or not this type can be 
observed in a water tank is left for a future work. 

We found that in most cases, although not always, 
a portion of the sand escaped from a transverse ripple 
before a barchan emerged. This sand escape can also be 
seen in the water tank experiments, although it is not 
clear whether or not emergence of a barchan is always 
accompanied by the sand escape, because observing it 
is difficult unless the amount is large enough. The sand 
escape is a precursor of a barchan formation, so that it 
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Fig. 9. Type B transition from a transverse ripple to a barchan. 
The barchan appeared in the windward side from fragments of 
transverse ripple. Initial condition for sand bed are Dq = 3 and 
Lo = 64. Time steps are (a) 65000, (b) 70000, (c) 75000 and (d) 
80000. 



Fig. 7. Dependence of ripple size on the length Lq. The depth 
and time steps were fixed as Do = 3 and t = 10000, respectively, 
(a) Lo = 128, (b) Lq = 512 and (c) Lo = 768. 
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Fig. 8. Type A transition from a transverse ripple to a barchan. 
The barchan appeared in the leeward side from fragments of 
transverse ripple. Initial condition for sand bed are Dq = 2.5 
and Lo = 256. Time steps are (a) 45000, (b) 50000, (c) 55000 
and (d) 65000. 



will be useful for locating the place where a barchan is 
emerging in natural dune fields. 

In conclusion, we succeeded in simulating a dune field 
dynamics by using a simple model in a steady unidi- 
rectional fiow. The model took only two processes into 
account: saltation and avalanche. The barchans emerged 
through transverse ripples from the initial sand bed. Be- 
fore barchans formed from the transverse ripples, a por- 
tion of sand often escaped from the ripples. 
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